Higher Moments of Net-Baryon Distribution as Probes of QCD Critical Point 
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It is crucially important to find an observable which is independent on the acceptance and late 
collision process, in order to search for the possible Critical Point predicted by QCD. By utilizing 
A Multi-Phase Transport (AMPT) model and Ultra Relativistic Quantum Molecular Dynamics 
(UrQMD) model, we study the centrality and evolution time dependence of higher moments of net- 
baryon distribution in Au + Au collisions at ^/sjvjv = 17.3 GeV. The results suggest that Kurtosis 
and Skewness are less sensitive to the acceptance effect and late collision process. Thus, they should 
be good observables providing the information of the early stage of heavy ion collision. In addition, 
our study shows that the Kurtosis times a 2 of net-proton distribution are quite different to that of 
net-baryon when collisions energy is lower than y/s^N = 20 GeV, the Monte Calor calculations on 
Kurtosis-o- 2 have a deviation from the theoretical predictions. 

PACS numbers: 25.75.Gz, 21.65.Qr, 24.60.Ky 
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I. INTRODUCTION 

In QCD phase diagram, the first order phase transition 
boundary separates the two phases: the hadron gas and 
the quark gluon plasma(QGP). The Critical Point(CP) 
which locates at the end of the transition boundary, is 
a distinct singular feature of QCD phase diagram[l|. To 
map the components of the QCD phase diagram is one 
of the main goals of heavy-ion physics, and searching for 
the critical point has been addressed in both theoretical 
and experimental studies 0-0], all these works help us to 
gather information of the singularities near the CP. If the 
CP existence around /is ^300MeV[|J it can be studied 
by heavy ion experiments, like BNL Relativistic Heavy 
Ion Collider (RHIC) with its Beam Energy Scan (BES) 
programJH, CERN Super Proton Synchrotron (SPS)@ 
and the future GSI Facility for Antiproton and Ion Re- 
search (FAIR)Q. 

The fluctuations of conserved quantities, like net 
baryon number, electric charge and strangeness, are con- 
sidered to be sensitive indicators for the structure of cre- 
ated system [§| . The characteristic feature of a CP is the 
diverge of the fluctuations. Most proposed fluctuation of 
observables are variations of 2 rd order moments of the dis- 
tribution, such as particle ratio [|| , charged dynamical @ 
and ^-measure It has been proved that all these 

observables are proportional to approximately £ 2 
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where £ is the correlation length which will diverges at 
the CP. Theoretical calculations predict that £ rj 2 — 3 fin 
for heavy ion collision [12], thus, it is extremely difficult 
to measure in experiments. However, the recent results 
shown that higher moments of conserved quantities dis- 
tribution are more sensitive to CP due to their strong 
dependence on £, the 3 rd order moment Skewness^ £ 4 - 5 
while the 4 th order moment Kurtosis~ £ 7 [||. Also the 



product Kurtosis- a 2 (called K eS in reference [l3j]) which 
is equal to the ratio of 4 th order to 2 nd order susceptibil- 
ities shows a large deviation from unity near the CP by 
lattice calculations and QCD based model studies [T3 |. 

Studying the event-by-event fluctuations [HI, [l|| is an 
effective way to address fluctuations of a system created 
in a heavy ion collision [l7j . In this paper, we will present 
the study of higher moments as a function of evolution 
time, and we will discuss the acceptance effect on higher 
moments of net baryon distribution. 



II. MONTE CARLO MODELS 

Monte Carlo event generators, AMPT and UrQMD, 
have been used in this study. A Multi-Phase Trans- 
port(AMPT) model is made up by four main parts: 
the initial conditions, partonic interactions, hadroniza- 
tion and hadronic rescattering. The initial conditions, 
which include the spatial and momentum distributions of 
minijet partons and soft string excitations, are obtained 
from the Heavy Ion Jet Interaction Generator (HIJING) 
model [l8j|. Scatterings between par tons are modeled 
by Zhang's Parton Cascade(ZPC) [19j . which presently 
includes only two-body scatterings with cross sections 
obtained from the pQCD with screening masses. In 
the default AMPT model(abbr. "AMPT Default") [H, 
partons are recombincd with their parent strings when 
they stop interacting, and the resulting strings are con- 
verted to hadrons using the Lund string fragmentation 
model [H [Hj]. In the AMPT model with string melt- 
ing (abbr. "AMPT StringMelting") [H, The transi- 
tion from the partonic matter to the hadronic matter 
is achieved by a simple coalescence model, which com- 
bines two quarks into mesons and three quarks into 
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baryons 24[ . The authors of AMPT model use hadronic 
cascade, which is based on the A Relativistic Trans- 
port(ART) model [25| . to describe the dynamics of the 
subsequent hadronic matter. Final results from the 
AMPT model are obtained after hadronic interactions 
are terminated at a cutoff time t cu t- When setting 
the time-step(in fm/c) for hadron cascade to 0.2 (de- 
fault value), the termination time of hadron cascade 
t cut = 0.2 x NTMAX. Here NTMAX is the number of 
time-steps, the default setting is NTMAX. When the 
time-step is fixed, larger NTMAX means longer time of 
hadron rescatterings. Note NTMAX lets all resonance 
rapidly decay(less than 0.6fm/c) and turns off hadronic 
interactions effectively In this paper, we measure the 
fluctuation which belongs to different processes: scenario 

(a) "parton" is a process between ZPC and quark coa- 
lescence, it still locates in the partonic phase, scenario 

(b) "w/o ART" is the time that hadronization(or quark 
coalescence) has finished, it is in hadronic phase but no 
hadron cascade happens, scenario (c)and (d)are origi- 
nated from different termination time t cut , by setting dif- 
ferent NTMAX value we can control the hadronic process 
in the simulated heavy ion collisions . It helps us to un- 
derstand the hadronic effect on the observable. In order 
to investigate different effects in time evolution, we will 
study the fluctuation of higher moments in each process. 
Comparing the measured value from different processes, 
we can study whether the fluctuation of the higher mo- 
ment can provide the information from the early stage of 
the collisions. 

The Ultra Relativistic Quantum Molecular Dynam- 
ics (UrQMD) model [2(| is also been used in this pa- 
per. It is a microscopic transport theory based on the 
covariant propagation of all hadrons on classical trajec- 
tories in combination with stochastic binary scatterings, 
color string formation and resonance decay. The UrQMD 
model represents a Monte Carlo solution of a large set of 
coupled partial intcgro-differcntial equations for the time 
evolution of the various phase space densities fi(x,p) of 
different species of particles. In the input file, one can 
control the time to propagate and output time-interval 
(in fm/c), in this paper we set the time to propagate 
is 40 fm/c and the time- interval is 2 fm/c. More detail 
descriptions can be found in Ref. (2r| . 
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FIG. 1: (Color online) The centrality dependence of Kurtosis 
of net-proton, net-cBaryon, net-Baryon distribution in differ- 
ent transverse momentum windows at ^/snn = 17.3 GeV by 
AMPT StringMelting. 
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FIG. 2: (Color online) The time dependence of Kurtosis of 
net-proton, net-cBaryon, net-Baryon distribution with b=2, 
\y\ < 0.5 and < p T < 1 (GeV/c) by UrQMD model. 



III. RESULTS AND DISCUSSION 

In experiments, neutrons can not be detected, and 
the reconstruction efficiency is relatively low for strange 
hadrons, especially for multi-strange hadrons, such as 3 
and Q [30]. Fortunately, theoretical calculation suggests 
that the proton number is a meaningful observable for 
the purpose of detecting the CP in heavy ion experi- 
ments [23|, its fluctuation completely reflect the singu- 
larity of the baryon number susceptibility. Thus, only if 
the measurements from net baryon and net proton distri- 
bution are similar, we can searching for the CP by mea- 



suring the fluctuation of various moments of net proton 
distribution. 

Observables that can reflect real dynamics and have 
little influence on the finite acceptance, are worth ex- 
ploring from the theoretical point of view [H, Ull • Due 
to the limited acceptances in experiments, it is necessary 
to study the acceptance effect by Monte Calor model. In 
the previous work |3lj . we presented that the Skewness 
and Kurtosis have little influence on the size of rapidity 
windows and transverse momentum regions, they could 
be good candidate for searching the CP in experiments. 

Firstly, we study the transverse momentum window 
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FIG. 3: (Color online) The centrality dependence of Kurtosis of net-proton, net-cBaryon and net-baryon distribution in time 
evolution at ^/Inn = 17.3GeV within AMPT StringMelting. 



cut effect on higher moments. Figure 1 shows the central- 
ity dependence of Kurtosis of net-proton, net-cBaryon, 
net-Baryon distribution at y / s]v7V = 17.3 GeV by AMPT 
StringMelting. The transverse momentum window are 
0< p T <0.4 (GeV/c), 0.4< p T <0.8 (GcV/c) and 0.8< 
Pt <1.2 (GeV/c), respectively. The only difference be- 
tween net-Baryon and net-cBaryon is that net-cBaryon 
excludes neutrons. As we presented before 3l|, the Kur- 
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tosis has a decreasing trend from peripheral to central 
collision. The value of Kurtosis approaches to zero in 
the central collisions, it means that these distributions 
are more similar to Gaussian distribution. Based on dis- 
cussions in Ref. [HJ, Kurtosis can be scaled by (A part ), 
if Kurtosis deviate from this (iV part ) scaling curve and 
exhibit a non-monotonic behavior, it would indicate the 
new physics which is related to CP. 

As shown in each panel of Fig. 1, the results from Net- 
Proton, Net-cBaryon and Net-Baryon are almost consis- 
tent with each other, especially in the central collisions. 
At the same time, the trends of Kurtosis in three trans- 
verse momentum windows seems similar, it suggests Kur- 
tosis has little dependence on the chosen transverse mo- 
mentum window. Together with the previous results of 
windows size dependence 3l|, we argue that Kurtosis is 



independent on acceptance. We also calculate the trans- 
verse momentum window cut effect on Skcwness and find 
that Skcwness is an observable which does not depend on 
acceptance. 

As an ideal observable for the signatures of CP, it 
should be only sensitive to the early stage of the colli- 
sions. Therefore, it is important to investigate the be- 
havior of these observablcs as a function of evolution 
time. In Fig. 2 we study the Kurtosis of net-proton, 
net-cBaryon, net-Baryon distributions as a function of 
evolution time by UrQMD model. The measured values 
of Kurtosis of net-Baryon, net-cBaryon and net-proton 
distribution keep a constant. The results suggest that 
the Kurtosis of the net-baryon distribution in the final 
state can reflect the distribution of the early stage. Most 
importantly, The system's evolution docs not affect the 
signature. 

Based on the UrQMD model calculation, we find the 
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FIG. 4: (Color online) Different particles yields in time evo- 
lution within AMPT StringMelting at ^/s^"=17.3 GeV. 



Kurtosis can provide the information of the signature of 
the early stage. With AMPT StringMelting, we can il- 
lustrate different effects on Kurtosis in time evolution 
clearly. In Fig. 3 we study the Kurtosis as a func- 
tion of centralities in scenario (a)parton, (b)w/o ART 
, (c)NTMAX = 3, (d)NTMAX = 150 with < p T < 
l(GeV/c) and \y\ < 0.5 at ^/snn =17.3 GeV, respec- 
tively. From scenario (a) to (b), the system experiences 
the hadronization process. Comparing panel (a) and (b) 
of Fig. 3, we find the Kurtosis of the distributions origi- 
nate from early partonic phase, the quark coalescence- 
like hadronization almost doesn't affect Kurtosis. In 
AMPT StringMelting, the hadron cascade procedure is 
regarded as two processes, resonance decay and hadron 
rcscattering. From scenario (b) w/o ART to scenario (c) 
NTMAX = 3, the system experiences the resonance de- 
cay process. Fig. 4 is the particles' yield in different 
evolution time in central collision at y / ijv r /v = 17.3 GeV 
by AMPT StringMelting. We can find all resonances de- 
cay before t cut = 0.6fm/c and the certain particle yield 
doesn't change after NTMAX = 3. Here the yield is the 
ratio of a certain particle number to the total particles 
number. When we compare panel (b) and (c) of Fig. 3, 
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FIG. 5: (Color online) The beam energy dependence of 
Kurtosis-cr 2 of net-proton, net-cBaryon, net-Baryon distribu- 
tions by AMPT StringMelting. 



it seems that the resonance decay process doesn't affect 
Kurtosis. 

On the other hand, large hadron rescattering effect 
may destroy the fluctuation which originates from the 
early stage of heavy ion collision. This effect depends 
on two factors: one is the time that particles go through 
the collision region, the other is the density in the colli- 
sion region [33j . In AMPT StringMelting, we can study 
the hadron rescattering effect on higher moments by con- 
trolling the termination time of hadron cascade(i cu t). 
NTMAX = 3 is regarded as that there is no hadron 
rescattering, while NTAMX = 150 is corresponding to 
termination time of 30 fm/c. From panel (c) and (d) of 
Fig. 3, it suggests the observed value of Kurtosis nearly 
does not change and almost keeps the same trend, even 
though the system experiences full hadronic rescatter- 
ings. It is reasonably to conclude that hadron rescatter- 
ings do not have clear influence on Kurtosis. 

Besides Kurtosis and Skewness, the product of 
Kurtosis and Variance (Kurtosis-cr 2 ) of the net-baryon 
distribution has also been used to search for the CP. This 
observable is related to the ratio of quartic (4 th order) 
and quadratic (2 nd order) cumulates of baryon number 
susceptibilities. Lattice QCD and some QCD based 
model [HI predicted that the fluctuation of Kurtosiscr 2 
changed rapidly at the transition temperature Tc- 
Furthermore, it is predicted that baryon number suscep- 
tibilities will diverge when close to the CP, this will bring 
about the deviation of Kurtosis-cr 2 from being constant. 
All of the works argue that Kurtosis -a 2 is a worthful 
observable. However, as we presented before [3l|, the 
Kurtosis-cr 2 is dependent on the acceptance, namely, 



lager acceptance lead to smaller measured value. Also, 
as shown in Fig. 5, results of Kurtosis-cr 2 of net-baryon 
distribution from two versions of AMPT model have a 
strong dependence on the collision energy, it is quite 
different from net-proton distribution. This difference 
may due to baryon stopping effect in low energy region. 
The similar result can be found in Ref. [l3[ by UrQMD 
calculation. 



IV. SUMMARY AND OUTLOOK 

In this paper we study fluctuation of higher moments of 
net-baryon distriubtion with AMPT and UrQMD mod- 
els. Based on these results, we find the higher moments 
of net-baryon distribution are independent on the cho- 
sen transverse momentum window. Together with the 
previous results of rapidity and transverse momentum 
windows size dependence, we conclude that Kurtosis and 
Skewness are independent on acceptance. 

We also study the time evolution effects and late 
hadronic effects on higher moments of net-baryon distri- 
bution in Au+Au collisions at -y/ijvlv = 17.3 GeV with 
UrQMD and AMPT model. The results seem that quark 
coalescence hadronization process, resonances decay pro- 
cess and hadronic rescattcrings don't affect the trend and 
value of Kurtosis and Skewness as a function of centrali- 
ties. Thus, they should be good oberservables to search- 
ing for the possible critical point predicted by QCD, non 
monotonic behaviors of the Kurtosis and Skewness as a 
function of beam energy or centrality will demonstrate 
the existence of a CP. The Kurtosis-cr 2 of net-proton dis- 
tribution show great difference to that of net-baryon dis- 
tribution when collisions energy is lower than ^snn = 
20 GeV, it deviates from the theoretical predictions. 

In the future, the RHIC Beam Energy Scan and GSI 
Facility for Antiproton and Ion Research will provide 
the possibility to locate the critical point in experiment. 
Higher moments of net-baryon distribution should be one 
of the powerful tools. 
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